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We show that the observed zenith angle dependence of 
the atmospheric neutrinos can be accounted for by neutrino 
decay. Furthermore, it is possible to account for all neutrino 
anomalies with just three flavors. 



I. NEUTRINO DECAY PHENOMENOLOGY 

According to the analysis by the Super-Kamiokande 
(SK) Collaboration 0| of their atmospheric neutrino 
data, the L/E dependence of the data is well ac- 
counted for by oscillation into (or a sterile neutrino 
Vst) with a mixing angle sin^ 29 ^ 0.8 and a Sm^ in the 
range 3 x 10~'*-8 x 10"'^ eV^. The i'^ survival probabihty 
in vacuum is given by 



(1) 



According to the data and in this simple two flavor 
mixing hypothesis, the i>e channel is unaffected and hence 
the survival probability Pee = 1- The interpretation of 
the atmospheric neutrino data in terms of oscillations is 
an old and venerable proposal 

It is obviously an important question to ask whether 
oscillation is the only possible explanation for the ob- 
served L/E dependence. Several other interpretations 
have been offered In this Letter we raise another 
possibility - neutrino decay. 

We begin by noting that decay implies a non-zero mass 
difference between two neutrino states and thus, in gen- 
eral, mixing as well. For definiteness, let us assume the 
existence of just three light neutrinos, and label as vi, 1^2, 
and i>3 that mass eigenstate with the largest admixture 
in the flavor state Ve, J^p, and Vri respectively. We further 
assume the dominant component of v^, i.e. 1^2, to be the 
only unstable state, with a rest-frame lifetime tq. There 
are strong limits coming from the nonobservation of tt 
and K decay to anomalous final states containing e^'s 
on the participation of i^e in non-SM vertices [Q. Con- 
sequently, i^e must nearly decouple from the unstable 1^2 
and its decay partner , and we are led to i^e vi, and 



cos l>2 



with 7712 > TT^3- From Eq. 
survival probability is 



mOiy^, (2) 
with an unstable 7^2, the 



P^^ sin"* d + cos^ 9 exp(-aL/S) 

+ 2 sin^ 9 cos^ 9 exp{~aL/2E) cos{Sm^L/2E) , (3) 



where Sm^ — 



m\ and a — 7712 /'Ha- 



lf, as we argue later, (5m^ > O.I eV^, then 
cofi{6m? L / 2E) effectively averages to zero for atmo- 
spheric neutrinos and P^^ becomes 



P^^ = sin'' 9 + cos^ 9 exp{-aL/E) . 



(4) 
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FIG. 1. The Super-Kamiokande data/expectations as a 
function of L/E, for electron events (upper) and muon events 
(lower). Our model normalized to the electron flux total is 
shown by the lines, indicating an acceptable fit for decaying 
muon neutrinos. 

We find that it is possible to choose 9 and a to provide 
a good fit to the Super-K L/E distributions of i/^ events 
and Vfji/ve event ratio. The fit to the L/E distribution 
of the events is shown in Figure |l|, where smearing 
over L and E is included. The best fit values of the 
two parameters are cos^ 9 ^ 0.97 (i.e. 9 ~ 10°) and 
a ^ 1 GeV/DE, where De = 12,800 km is the diameter 
of the Earth. This best-fit a value corresponds to a rest- 
frame 7^2 lifetime of 
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To = m2/a 



1712 

(1 eV) 



X 10 



-10 



(5) 



Such a lifetime and decay length are not in conflict 
with current limits on non-radiative modes . Previous 
limits on decay lengths from accelerators are only of 
the order of a few km. 



II. NON-RADIATIVE NEUTRINO-DECAY 
MODELS 

The simplest possibility for fast invisible neutrino de- 
cay is to choose 1/2 and to be Dirac states, with a de- 
cay interaction between the right-handed 5f7(2)-singlet 
states given by Q 



jOi = 523 I'sfl i^2fl X + h.c. 



(6) 



where x is a complex scalar field with lepton number 
—2, Iw = 0, and hypercharge zero, x is not a Nambu- 
Goldstone boson. The scalar field should be light com- 
pared to the neutrino masses. 

The quantum numbers of x make it difficult to test its 
existence in laboratory experiments. Processes such as 
H —>■ e + X are forbidden. Furthermore, neutrino decay 
is not helicity-suppressed but the interactions of its decay 
products are: all new processes involving neutrinos are 
suppressed by the chirality factor (m^/E^)'^. 

A decay model for Majorana neutrinos can also be con- 
structed. The effective decay interaction is 



-C/ = 523 J + h.c. . 



(7) 



where J is a Majoron field which has to be dom- 
inantly singlet with only a small triplet admixture, in 
order to satisfy the constraints from the invisible decay 
width of Z. A full model needs futher embellishment to 
generate the above coupling at tree level and generate 
the desired neutrino masses and mixings. The model de- 
scribed by Acker et al. can be easily adapted to our 
purpose here. 

With the interaction of Eq. or Eq. (|^), the rest- 
frame lifetime of 1^2 is given by 



To 



IGtt 



5I3 5m? {m2 + m^Y 



(8) 



and hence 

5I3 5m? — IGttq! (1 -f 771,3/7712)^^ , 
leading to, for < 7773/7712 < 1, 

5I3 W - (2-7) X lO^'^eV^. (10) 

From studying K decays, in particular the modes K 
11 + neutrals, a bound on the coupling 523 can be derived 
which is IHl 



< 2.4x10-4, (11) 
thus leading to a bound on 5m'^ in the present model 



of 



5ml^ > 0.73 eV^ . 



(12) 



This result justifies the above approximation of large 
5m'^. 

In both the Dirac neutrino and the Majorana neu- 
trino decay models, the small mass scale for the scalar 
fields seems to require fine-tuning. This is unavoidable, 
though unappealing. The Dirac and Majorana models 
are testably different in that in the Dirac-neutrino model 
the decay is into essentially undetectable final states, 
whereas the Majorana-neutrino model posits decay of 
I'^'s and 7^'s into i^r's and v^-^s respectively, which are ob- 
servable in principle. Furthermore, with the Majorana- 
neutrino model, lepton number is broken and neutri- 
noless double-beta decay is allowed; with the Dirac- 
neutrino model, lepton number is unbroken and neutri- 
noless double-beta decay is not allowed. 

If the decay of the V2 is into a new(sterile) neutrino 
with which it does not mix, then the 5rn? appearing in 
the oscillation is not restricted by the constraint of Eq. 
(12) and can be very small. Then the cosine term in 
Eq. (3) is essentially unity and the survival probability 
is given by 



P^^ = [sin^ e + cos^ 6iexp(-aL/2£;)]^ 



(13) 



In this case, even better fits to the Super-K data can 
be obtained. However, we do not pursue this class of 
models further at this time. 



III. INCLUSION OF ALL NEUTRINO 
ANOMALIES 



It is easy to incorporate the solar and LSND neutrino 
anomalies into the present discussion. The decay of V2 
with Sm^'^ ^ 0.73eV^ explains the atmospheric data, but 
also allows 5m\^ to accommodate the LSND result. We 
are free to choose the remaining 5w? to accommodate the 
solar anomaly. As mentioned earlier, Ues must be rather 
small; and hence the solution to the solar v^. depletion is 
necessarily small-angle mixing enhanced by the MSW ef- 
fect. To this end, we set Sm?,, 



5ml^ ~ 0.5 X 10"^ eV^ 



to complement the (Jmjs of Eq. (|12| 



Atmospheric v 's 



With these 5mjj, the Vfj. survival probability is given 
by Eq. M) with cos^ 9 = \U^2\^ - 0.97 and sii? = 
1 — I^^mzI 7 ^ii'i the earlier two-flavor discussion and flt 
remain intact. The survival probability for v^s is 



Pee=(l-|C/e2n' + |C/e2pe-"^/^ 



2|C/e2|'(l - \Ue2f) cos{5ml2L/2E)e 



-ahjlE 



(14) 
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If Ue2 ^ 1, then Pee ~ 1, as observed. 
Solar V 's 

For solar L/E values, the 1/2^8 have decayed away and 
the Ve survival probability is given by 



P.. 



1 



If/, 



iC/esl'sin^ 



^31 



6m\iL 



(15) 



By choosing 4|C/ei|^|J7e3p w 5.5 x lO^'^, with TO3 > 
mi, one can reproduce the small-angle MSW solution for 
solar neutrinos H. Furthermore, the resulting value of 



U, 



e3l 



1.4 X 10 easily satisfies the upper bound on the 
coupling which is gl = gIziW e?,]"^ + |t^e2p) < 3 x 10~^ 

i- 

LSND 

At the L/E value relevant to the LSND experiment, 
the v^-Ve conversion probability is given by 



P^e =4|C/,,2r|f/e2rsin- 



AE 



(16) 



With (5m^3 ~ 0(1 eV)2 and |J7^2p - 1, choosing 
p 

/ie 



|C/e2p ~ 10 ^ allows this P,jk to account for the LSND 



observations 

Summary of Mass and Mixing 

To summarize, the three-neutrino mixing matrix with 
the approximate form 



0.999 0.02 -0.04 
\ -0.05 0.985 0.17 
e -0.17 0.985 



(17) 



and mass differences given by 

(5m^i « (5m^3 > 0.73 eV^ , and 5m?3 - 0.5 x 10"^ eV^ , 

(18) 

along with an unstable V2 satisfying the lifetime con- 
straint of Eq. (|^), explain all three neutrino anomalies 
without violating any known data. The very small entry, 
e in Eq. (p7|), is fixed by unitarity of the mixing matrix. 



Big Bang Nucleosynthesis 

A neutrino with lifetime as short as the one dis- 
cussed here decays before the primordial neutrinos de- 
couple (about 2 MeV). Then the decay products would 
also achieve thermal equilibrium, increasing the effective 
number of neutrino species to nearly 5, and would 
be in conflict with the bound on light degrees of freedom 
which is generally considered to be about 3.5 |12| or even 
lower ]l3t . Given all the uncertainties, more cautious and 
conservative estimates |^4j lead to upper bounds on N^, 
in the range of 4 to 5. Furthermore, the bound is mod- 
ified if there is lepton number asymmetry present p5| . 
Hence, although this may be a potential problem, we do 
not regard it as fatal. 



Supernova Emission 

The decay of at the fast rate envisaged here will cer- 
tainly modify supernova dynamics, since the decay will 
occur inside the core. In the Dirac case, the sterile de- 
cay products will carry away energy on a very short time 
scale and this probably conflicts with the observed period 
of a few seconds of the SN1987A neutrino burst. In the 
Majorana case, the decay products are not sterile and 
this problem does not arise, and the effects on dynamics 
should be much milder. It is not clear what is "the pre- 
cise range of parameters that can be ruled out or ruled 
in by the SN1987A signal" |l6|. 

Cosmic V Fluxes 

In neutrinos coming from distant sources, such as Su- 
pernovae, AGN's, and GRB's, the v^{~ i^2)'s have de- 
cayed away and only Vf,,Vf,^Vr and Vt will arrive at the 
Earth. Existing neutrino telescopes, as well as those un- 
der construction, will not observe the tiny component 
(^ sin** 10^'') of surviving j/^'s (apart from atmo- 
spheric ones). 

Future Tests 

There are several opportunities to test the neutrino 
decay hypothesis decisively. One is that the Super- 
Kamiokande collaboration, with sufficient data, can dis- 
tinguish between the two LjE distributions: the oscil- 
latory one as given in Eq. (|l^) and the decaying one as 
given in Eq. (^j. A second is that in the upcoming /i's in 
Super-K (which come from i^^'s of average energy around 
100 GeV) the decay should have very little effect due to 
the long lab-frame lifetime of the 1^2 and the zenith an- 
gle distribution of events should be undistorted. For the 
same reason, partially contained /i's should show some 
up-down asymmetry but less pronounced than for the 
fully-contained events. 

A third opportunity exists with future long-baseline 
experiments |17|| . Expectations are quite different with 
the neutrino decay interpretation of Super-K results, 
compared to the oscillation interpretation. In the de- 
cay scenario the typical v^^ survival probability is 94%, 
the v^-Vt conversion probability is about 6%, and the 
v^j^-Ve conversion probability about 0.2%. (There is, in the 
Majorana model, a further small Vt flux from the decay 
products; however their energies are much lower than the 
parent I'^'s due to the backward-peaked nature of the de- 
cay.) 

Conclusion 

We have presented a neutrino decay scenario capable of 
explaining the atmospheric neutrino data. Furthermore, 
it enables us to explain all neutrino anomalies with just 
three neutrino flavors. The scenario violates no available 
data, and is economical in that only a single physics pa- 
rameter {a — 1712 /tq) is added beyond those of the usual 
oscillation phenomenology. The decay possibility can be 
checked in future Super-K data as well as in forthcoming 
long-baseline experiments and neutrino telescopes. 
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